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Corporate defaults may be triggered by some major market news or events such as fi- 
nancial crises or collapses of major banks or financial institutions. With a view to develop 
a more realistic model for credit risk analysis, we introduce a new type of reduced-form 
intensity-based model that can incorporate the impacts of both observable "trigger" events 
and economic environment on corporate defaults. The key idea of the model is to aug- 
ment a Cox process with trigger events. Both single-default and multiple-default cases 
are considered in this paper. In the former case, a simple expression for the distribution 
of the default time is obtained. Applications of the proposed model to price defaultable 
bonds and multi-name Credit Default Swaps (CDSs) are provided. 
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1 Introduction 



Modeling default risk has long been an important problem in both theory and practice 
of banking and finance. In the aftermath of the global financial crisis (GFC), much 
attention has been paid to investigating the appropriateness of the current practice of 
default risk modeling in banking, finance and insurance industries. Popular credit risk 
models currently used in the industries have their origins in two major classes of models. 
The first class of models was pioneered by Black and Scholes (1973) and Merton (1974) 
and is called a structural firm value model. The basic idea of the model is to describe 
explicitly the relationship between the asset value of a firm and the default of the firm. 
More specifically, the default of the firm is triggered by the event that the asset value 
of the firm falls below a certain threshold level related to the liabilities of the firm. The 
structural firm value model provides the theoretical basis for the commercial KMV model 
which has been widely used for default risk model in the financial industry. The second 
class of models was developed by Jarrow and TurnbuU (1995) and Madan and Unal (1998) 
and is called a reduced-form, intensity-based credit risk model. The basic idea of the model 
is to consider defaults as exogenous events and to model their occurrences using Poisson 
processes and their variants. In this paper, we focus on reduced-form, intensity-based 
credit risk models. 

Reduced-form, intensity-based credit risk models have been widely used to model 
portfolio credit risk and to describe dependent default risks. There are two major types 
of reduced-form, intensity-based models for describing dependent default risk, namely 
bottom-up models and top-down models. Bottom-up models focus on modeling default 
intensities of individual reference entities and their aggregation to form a portfolio default 
intensity. Some works on bottom-up models include Duffie and Garleanu (2001), Jarrow 
and Yu (2001), Sch5nbucher and Schubert (2001), Giesecke and Goldberg (2004), Duffie, 
Saita and Wang (2006) and Yu (2007) etc. These works differ mainly in their specifications 
for the parametric forms of default intensities of individual entities and the way these 
intensities are aggregated. The top-down models concern modeling the occurence defaults 
at a portfolio level. A default intensity for the whole portfolio is modeled without reference 
to the identities of individual entities. Some procedures such as random thinning can be 
used to recover the default intensities of the individual entities. Some works on top- 
down models include Davis and Lo (2001), Giesecke, Goldberg and Ding(2011), Brigo, 
Pallavicini and Torresetti (2006), Longstaff and Rajan (2008) and Cont and Minca (2011). 

We focus on the bottom-up model. Lando (1988) proposed a reduced-form, intensity- 
based model, where the occurrence of a default is described by the first jump of a Cox 
process. The main advantage of the Lando's model is that under his model, a simple 
pricing formula for a defaultable risky asset can be obtained. This formula is similar 
to the one for the default-free counterpart of the risky asset. Yu (2007) extended the 
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Lando's model to incorporate multiple defaults and their correlation. The so-called "total 
hazard construction" by Norros (1986) and Shaked and Shathanthikumar (1987) was used 
to generate default times with interacting intensities. Zheng & Jiang (2009) proposed a 
unified factor-contagion model for modeling correlated defaults and provide an analytical 
solution for modeling default times with "total hazard construction". Gu et al. (2011) 
introduced an "ordered default rate" method to give a recursive formula for the distri- 
bution of default times in pricing basket Credit Default Swaps (CDSs) in the context of 
a reduced-form, intensity-based model, which significantly enhances the computational 
efficiency in finding the prices of CDSs. One of the shortcomings in a number of exist- 
ing reduced-form intensity-based models is that they fail to incorporate the impact of 
major market events, such as financial crises, on corporate defaults. This may lead to 
underestimation of default risk and also undervaluation of defaultable risky products. 

In this paper, we address the problem of how to incorporate the direct impact of ob- 
servable "trigger" events on corporate defaults in a reduced-form, intensity-based credit 
risk model. The key idea is to describe the occurrence of a default by the first unrecover- 
able "trigger" event. Armed with a Cox process for default risk, we incorporate the impact 
of economic environment on defaults by allowing the default intensity depending on an 
underlying state process representing the variation of economic environment over time. 
We consider both the single-default and multiple-default cases. In the single-default case, 
we obtain a simple expression for the distribution of the default time. This distribution is 
useful for pricing defaultable securities. We then extend the model to the multiple-default 
case with a view to incorporating default correlation. To provide a tractable and practi- 
cal way to value defaultable securities, we focus on the case where the state process for 
economic environment is modeled by a continuous-time, finite-state, observable Markov 
chain. Applications of the proposed model to value defaultable bonds and basket Credit 
Default Swaps (CDSs) are discussed. We also provide numerical results to illustrate the 
sensitivity of the prices of these securities with respect to changes in key parameters. 

The rest of the paper is organized as follows. Section 2 presents the basic model 
and derives the distribution of a default time. Section 3 provides the extension of model 
framework to the case of multiple correlated defaults. Section 4 presents the Markov 
chain model for the state process of economic environment. Applications to pricing the 
defaultable securities are given in Section 5. We then conclude the paper in Section 6. 

2 The Basic Model 

A popular reduced-form, intensity-based credit risk model was proposed by Lando (1998), 
where the occurrence of a default was described by the first jump of a Cox process with 
stochastic intensities {Xt}t>Q depending on an underlying state process {Xt}t>o describing 
the evolution of economic environment over time. Here we aim at extending the Lando's 
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model by incorporating explicitly the "trigger" events such as financial crises or extra- 
ordinary market news into default risk modeling. We assume that these "trigger" events 
are observable and may lead to the default of a corporation. Furthermore, we suppose 
that the corporation may recover from a "trigger" event via re-organizing its resources or 
re-structuring. We assume that the emergence of "trigger" events are modeled by a Cox 
process, which is also called doubly Poisson process in the statistical literature and has a 
remarkable history in statistics. In what follows, we shall describe the mathematical set 
up of the Cox process describing the "trigger" events. 

Uncertainty is described by a complete filtered probability space {fl, J^, {J^t}t>o, P), 
where P is a given probability measure and {Tt}t>o is a filtration satisfying some usual 
conditions, (i.e., the right- continuity and the P-completeness). We shall define precisely 
the filtration {J^t}t>o in later part of this section. 

To describe the evolution of the state of economic environment over time, we define 
a state process {Xt}t>o- We assume that {Xt}t>o is a cadlag, (i.e., right continuous with 
left limits), process on (f2, J^, {J^t}t>o,P) with state space JJ. Let {A't}t>o be a standard 
Poisson process on (fi, J-", {J-'t}t>o, P), with Nq = 0, P-a.s. Write {Xt}t>o for a bounded, 
non-negative stochastic process on {Q, J^, {J^t}t>o, P)- We assume that for each t > 0, 
Xt := X{Xt), for some non-negative continuous function A and that {Nt}t>o and {Xt}t>o 
are stochastically independent under P. For each t > 0, we define the following cumulative 
process {At}t>o- 



For each z = 1, 2, ■ ■ -, let r* be the arrival time of the i*^^ "trigger" event, which is modeled 
as the arrival time of a jump in the Cox process. Once a "trigger" event occurs at time s, 
a loss occurs to the firm , which is modeled as an arbitrary independent random variable 
"L ". We write {Ct}t>o ioi the process depending on the state process {Xt}t>o, where 

^ When we wish to evaluate the risk of a credit portfoho, we need to use a real-world probability 
measure. In this case, P can be interpreted as a real-world probability measure. On the other hand, 
when we wish to price defaultable securities, we must use a risk-neutral probability measure. In this case, 
there arc two approaches to interpret P. The first approach is to interpret P as a risk-ncutral probability 
measure and start with the risk-neutral probability measure directly. The second approach is to interpret 
P as a real- world probability measure and then use a measure change for Poisson processes to transform 
the real-world probability measure to a risk-neutral one. To simplify our discussion, when we discuss the 
pricing of defaultable securities, we shall adopt the first approach. 

^ Strictly speaking, the intensity measure A is defined on the cr-field generated by bounded subsets of 
the interval [0,cxd). The intensity measure A also has a Random density, or Radon-Nikodym derivative, 
A; by definition. 




Then a Cox (point) process {Nt}t>o with intensity measure 



A := {At}t>o a is defined by: 



Nt := Na, . 
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Ct := C{Xt) as the threshold value at time t, for some non-negative continuous function 
C. If L < Csi then the firm can recover from the "trigger" event; otherwise, the firm 
defaults. Davis & Lo (2001) introduce a Bernoulli contagion variable in a homogeneous 
setup, which is similar to variable "L" Here. 

Let {r*}j=i^2, -- be a sequence of stopping times representing the arrival times of "trig- 
ger" events defined by 

t' ■= mi{t > : iVt > z} , 

and {U}i=i^2,- - a sequence of arbitrary independent and identically distributed random 
variables. We assume that {r*}j=i 2,... and {L*}j=i^2, -- are stochastically independent under 
P. Define a random variable K taking values in {1, 2, ■ ■ ■} by 

K := min{2 : U > C^i} . 

Then the default time of a firm r is defined by: r := . We now specify more explicitly 
the information structure of our model. Define the filtrations {Gt}t>o, {'Ht}t>o and {Xt}t>o 
as follow: 

Gt := a{Xs : < s <t}y Af, 
Ut := a{Ns : < s < t} VAT , 

and 

It := a{l{r<s} : < s < t} V A/", 

where M is the collection of all P-null subsets in J-" and 1^ is the indicator function of 
an event A. Here we assume that the filtration {J^t}t>o is specified as follows: J-i := 
Qt V "Hj VXj. This represents the full observable information structure in our model. 

For each s > 0, let be the probability that the firm can recover from a "trigger" 
event if the event occurs at time s given the underlying state Xs. Then, 



Let Ps '■= ^ — Is- Then 



:= P{L < C{X,) I X,) 



Ps = P{L > C{X,) I X, 



The following result is one of our main results which gives the conditional and uncondi- 
tional distributions of the default time r. The proof can be found the Appendix. 

Proposition 1 For any t > s > 0, 

P{t > s I ^i) = exp | - ^ puKdu^ 



and 

P{t > s) = E 
where E is an expectation under P. 



exp PuKdu 
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The following results are also important to characterize the probability laws of the 
default time r. Their proofs can be found in the Appendix. 

Lemma 1 For any s < t, 

E{l{r>t} I VHs V Is) = l{r>s} exp |- ^ PuKduj . 
Lemma 2 The process 

^{t>T} ~ Pu^vX{u<T}du , t > , 

is an {{J^t}t>i)i'P) -martingale. 

To price the defaultable securities in the proposed model as above, we construct three 
"building blocks" as in Lando(1988). Before that, we assume that all of the expectations 
in this paper are taken under an equivalent martingale measure. We suppose that T 
denotes the expiry date of all contingent claims. The three "building blocks" are as 
follows: 

(I) A payment X G at a fixed date T which occurs if there has been no 

default before time T. 

(II) Ysl{T>sy- A stream of payments at a rate specified by the {^f}f>o-adapted process 
Y which stops when default occurs. 

(Ill) Zt-: a recovery payment at the time of default, where Z is a {^t}t>o-adapted process. 

Now we proceed to give the pricing formula of these three "building blocks" . 

Proposition 2 Suppose 

exp{-/o^r,ds}X, 

/o^ Ys exp{- /o rudu}ds, 

Jq ZsXsPs exp{- /o (r^, + \uPu)du}ds 

are integrable random variables. Then, 

E{exp{- rJs}Xl{,^T}) = ^(exp{- /o^(r, + p,A,)rfs}X), (1) 

E{J^Ysl{r>s}exp{-jQrudu}ds) = E{Jq Y,exp{- J^{ru + PuK)du}ds) (2) 

and 

E{exp{- Jq rsds}Zr) = E{Jq Zs\sPsexp{- Jq^Tu + XuPu)du}ds). (3) 

Proof: The results follow directly from Proposition 1. n 
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3 An Extension to Dependent Multiple defaults 



In this section, we shall extend the basic model in the last section to the multiple-default 
case. Again we consider the filtered probability space J-", {J-i}t>o, -P), where the fil- 
tration satisfies the usual conditions and is specified as in the last section. Here the 
underlying state process of economic environment {Xt}t>o is also defined as in the last 
section, (i.e., a cadlag process). For each i = 1, 2, ■ ■ ■ , n, let Tj be a stopping time repre- 
senting the default time of the i^^ individual obliger. 

Let {N]:}t>o, i = 1,2, ■ ■ ■ ,n,he n independent standard Poisson processes with Nq = 0, 
P-.a.s. Write, for each i = 1,2, {Xl}t>o for a bounded, non-negative stochastic 

process which is adapted to the enlarged filtration containing the filtration {Qt}t>o 
the filtration generated by the sequence of stopping times {rj}i=i,2,...,n, where the latter 
filtration will be defined precisely later in this section. 

For each i = 1,2, ■ ■ ■ ,n and each t > 0, we define the intensity measure A* := {Al}t>o 
for the z^'^-individual obliger by putting: 

A! := / Vds < oo . 
Jo 

Then for each i = 1,2, , ■ ■ ■ ,n, a Cox (point) process {Nl}t>o for the i^^ obligor asso- 
ciated with the intensity measure Aj is defined by: 

:= iVlj . 

Define the arrival time of j^^ "trigger" event of name i, 

ri =M{t:Ni>j}, 1 = 1,2,. ..,n. 

Consider an array of i.i.d. arbitrary random variables 

{Ll,t = l,2,...,n,j = 1,2,...}. 

We suppose that this array of random variables is independent of all previous events. 
Define the random variables: 

Ki := mm{j : > C^,} , i = 1,2, - ■ ■ ,n, 

i 

where C* := C^{Xs) is the threshold value of name i at time s and a non-negative 
continuous function for z = 1, 2, . . . , n. 

Then the default time of name i is defined as 

'« • 'i • 

The information structure of the multiple-default model is specified as follows: 

nt ■■=a{iN:)ti:0<s<t}VAf , 
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and 



It:=a{il{r,<s})7=i-0<s<t}yAf. 



Again we assume that 



We also define 



pl := P{L > C\X,) I X,), ql : 



We further assume that for each i = 1,2, ■ ■ ■ ,n, P(0 < < oo) = 1 and P(rj = r,) = 
for any i ^ j. Moreover, X is an "exogenous" stochastic process in the sense that Qoo 
and Tit VX^ are conditionally independent given Qf 

We remark that this framework for multi-name defaults allows the default intensities 
to be sensitive to the observed defaults as well as the underlying state process X. Taking 
the consideration of the real practice that once a firm recovers from a "trigger" event, 
the impact of such a "trigger" is restricted to a very minor level, {\t}t>o is supposed 
{Qt VXt}t>o-adapted. 

For each i = 1, 2, ■ ■ ■ , n, if we treat information about the state process X and observed 
defaults of other names up to time t as a new "^f", then applying Lemma 2 in the single- 
name default case, we deduce that the process defined by: 



is an ({J-f }t>o, P)-martingale. 

Consequently, the total hazard construction method, Yu(2007) and Zhang & Jiang 
(2009), the order default rate approach by Gu et al.(2011) can be applied to compute 
the multi-name default time distribution under this framework. We also remark that this 
framework is indeed a generalization of the standard reduced-form intensity-based model 
by Lando (1988) and Yu (2007). If we set pl = 1, i.e., the firm cannot recover from 
"trigger" event P-a.s., then our proposed model is reduced to be the standard reduced- 
form intensity-based model. 

4 State Process as an Observable Markov Chain 

Suppose now that the state process {Xt)t>o follows a continuous-time, homogeneous, 
M-state Markov chain on {Q,J^,P) with state space {xi,X2, ■ ■ ■ ,xm}- The amount of 
time the chain X spends in state Xi before making a transition into another state is 
exponentially distributed with rate Vi . When the process leaves state Xi, it next enters 
state Xj with probability pij, where 




Pa = and "^Pij = 1 
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for all i. 

Let Tij{t) be the occupation time of the chain X in state Xj in the time interval [0,t] 
starting from Xo = Xi. We wish to determine the joint distribution of (Tji(t), Tj2(t), . . . ,TiM(t)). 
Note that the joint distribution of (Tji(t), Tj2(t), . . . ,TiMit)) is completely determined by 
its joint moment generating function. We shall derive the joint moment generating func- 
tion in the sequel. 

For each z = 1, 2, ■ ■ ■ , M, let 



T,(t) = (Ta(t),T,2(t),...,T,M(t))' 



and 



The moment generating function of Tj(t) is given by: 

^,{u,t) = E{exp{u^T,{t)}) 

Let denote the time of the first jump from state Xj to another state, i.e., ~ exp(fj 
Hence 

^iiu,t) = E{exp{u^T,{t)}) 

= E{E{exp{u^T,it)} I e.)) 

= T.P^kV^ f e("»-''»)^E(exp{iz^Tfc(t - s)})ds + e("»-''»)* 



k^i 

YpikVi f e("'-^')(*-^)^fc(M, s)ds + e^"'-''")*. 



k^i 

Taking the partial derivative with respect to t on both sides yields, 
d 



dt 



(4) 



For simplicity, we write 



= (^i(M,t),^2(M,t),...,^M(M,t))' 



and define the matrix 



A 



Ul - Vl P12V1 ■ ■ ■ PimVi 

P21V2 U2-V2 ■■■ P2MV2 



PmiVm PM2Vm ■ ■ ■ um - Vm 

Hence Equation (jl]) can be rewritten as a system of homogeneous linear differential equa- 
tions with constant coefficients 

^J{t) = A^^it), (5) 

where \E'm(0) = 1. Solving the ODEs, by using the Fundamental Theorem for Linear 
Systems (see Chapter 1 in L. Perko (2001)), yields the following proposition. 
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Proposition 3 The moment generating function of T.i(t) is given by \l/i(-u,t), where 
/ias a unique solution as 

This result will be used to price basket credit defaut swaps under our proposed model in 
Section 5.2. 

5 Applications 

In this section, we apply the ordered default rate approach in Gu et al. (2011) to price 
different defaultable securities under various assumptions of default correlation structures. 

5.1 Defaultable Bonds 

We first discuss the pricing of defaultable zero-coupon bonds with zero recovery under a 
constant default-free rate r in the case of two firms. The defaultable bond price is then 
proportional to the conditional survival probability P{Ti > T \ J^t). A "looping default" 
case with two firms was presented by Jarrow and Yu (2001), where they assumed that 
the default times of the two firms and have the following intensities: 

\f = ai + a2l{t>TB} 

and 

Af = 6i + 62l{t>T^}- 

Yu (2007) gave the default time distribution of this two-firm case in the standard reduced- 
form intensity-based modelling framework. In this section, we provide a solution of the 
two-firm case in our proposed modeling framework, where p\ = p. In other words, the 
firm recover from a "trigger" event with a constant probability 1 — p. 

Let A denote the first default time of these two names, the the default rate of 
-g _|_ gy Proposition 1, 

P(r^ Ar^ >t) = exp{-p{ai + bi)t). 

H6I1C6 

ai + hi 
and also 

F(r^ — > t \ > T^) = e~p("i+''2)t 
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Then by making use of convolution, we have 

dt ~ ' hi — a2 

dt ~ ' ai — 62 



Therefore, the marginal density of r and r are given by 



bi{ai + a2)p 
hi - 02 



ni — n.n ^ ' 



-p(ai+fei)t 



ai - O2 ^ ^ 
Comparing the results in Yu (2007), the marginal density of and under the standard 
reduced-form intensity-based model are, respectively, given by: 

9tA^) = fr^it, 1) and g^B{t) = f^sit, 1). 

Hence as the remark in Section 3, our proposed model is a generalization of the standard 
one. 



5.2 Basket Credit Default Swaps 

In this section we discuss the pricing of a basket CDS in the context of the multiple-default 
model described in Section 3 and provide the sensitivity analysis of the CDS premium by 
varying key model parameters. We consider a basket CDS contract that pays $1 if kih.- 
to-default out of a portfolio of reference entities occurs prior to expiry date. To simplify 
our discussion, we assume that the payment (if any) occurs at expiration, and that the 
buyer pays a premium at the initiation of the swap contract. With a constant risk- free 
interest rate r, the premium on the /cth-to-default CDS is given by: 

Sk = exp{-rT}P(r^' < T) , 

where is the /cth-to-default time and T is the expiry date. 

We assume the following default intensity for n reference names: 

Ai = X,(l + 6^1{,^.<i|) , ^ = l,2,---,n, 

and for all p\ = 1 — e"*^"^*. 

Let r'^ denote the fcth default time and is assigned to be 0. Let 

\^^^ =Xu{l + h{k-l)){n-{k-l)) 

denote the /cth-to-default rate as in Gu et al. (2011). By Proposition 1, we have the 
following relations: 

P (r'' - r'--i > s I a(r'--i) V = exp |- X^ubudu\ . 
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Hence we have 



Then by convolution, given Q^o, we have the following recursive formula for the conditional 
density function of r'^: 

frk\g^{t) = ^f^Pt /^fe-i|g^(u)exp|-^ X'^^^psds^du. 

Let Yt = XtPt and by the recursive formula, assuming b ^ l/i for alH = 1, 2, . . . , n — 1, 
the PDF and CDF of respectively, given the evolution of {Xt)t>o, are given by 

P{r' < t I G^) = ^(1 - exp{-/3, /* Y^du}) 
where the coefficients are given by the following recursion: 

Pk 

ak+i,j = \ 



(3k 

k-l 

- «fc+i,«, j = k 

u=0 



j = 0,l,...,A;-l 



P, = {n-j){l+jb) 



where ai o = n. Hence, 



k-l .1 
frkit) = J2'^k,jE{Ytexp{-l3j / Yudu}) 

n—n •'0 



3=0 

k—1 f 

Pir' < t) = E ^(1 - ^(exp{-/3, / Y^du})). 
j=o Pj •'^ 

If we assume {Xt)t>Q is a continuous time time-homogeneous Markov chain having M 

states {xi,X2, . . . ,xm}- By using the results in Section 4, we find the swap premium 

Sk- Indeed, Yt = XtPt is also a continuous time time-homogeneous Markov chain having 

M states {yi,y2, ■ ■ ■ yUAi}, where tji = Xi{l — e~^-^*). Then starting from state Xj, let 

y = (2/1, 2/2, • • • , yuY, we have 

P{r' <t) = E ^(1 - i?(exp{-/3,/T,(t)})) 
j=o Pj 

j=0 Pj 
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Consequently we have 

j=o Pj 

To give an example of the pricing of basket CDS, we set the number of states of (Xt)t>o, 
M = 4, the states Xj = O.lz, pij = 1/3 for any i ^ j, v = (wi, f2, ^3, ^^4) = (3, 2, 1, 3). We 
remark that the four states can represent the macroeconomic environment as "prosper- 
ous", "good", "neutral", "bad". We assume there are 10 firms in the industry {n = 10) 
and we start from state Xi, ie, Xq = xi, the default free interest rate r = 0.05, the expiry 
date T = 5 years. To examine the effect of contagion parameters b and c, we present the 
swap premium with the contagion parameters varying. 

As shown in Figure 1, an increase in the parameter b indeed raises the premium 
for all k as our intuition while the first-to-default swap price remains unchanged due to 
the contagion has no effect for the first-to-default time. The premium increases as the 
parameter c becomes bigger since the recovery from a "trigger" becomes more difficult in 
a weak macroeconomic environment. 

6 Concluding Remark 

In this paper, we introduce an extended reduced-form intensity-based model with "trig- 
ger" events, which captures an important feature in the market, namely the observable 
events that trigger defaults. Furthermore, we extend the model into a multi-name default 
model, with intensities driven by the history of the exogenous state process represent- 
ing the macroeconomic environment and the observed defaults. A Markov chain model 
for the state process is also proposed to model the macroeconomic environment and the 
distribution of the occupation time of each states is deduced by solving a system of ho- 
mogeneous ODEs. We demonstrate the pricing of two defaultable contingent claims using 
the ordered default rate approach by Gu et al. (2011) with numerical examples. 

There are still many outstanding issues for further research such as multi-state mi- 
gration model that links loss sizes of a firm due to "trigger" events and its credit states 
which in turn affect payoffs to investors, and hybrid model of "trigger" events and default 
contagion with effective computation. We are currently investigating these problems. 
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The basket CDS premiums with c=1 and varyiing b 




5 6 
kth to default 



b=0 

b=0.05 

b=0.1 
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The basket CDS premiums with b=0.05 and varying c 




4 5 6 7 
kth to default 
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Figure 1: The effect of the contagion parameters on the basket CDS premiums 
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7 Appendix 



7.1 Proof of Proposition 1 

Proof: Note that 

oo 



i=l 

oo 



1=1 

oo 

= X! -^[l{T«<s}PriQ'r»-l9r'-2 • • • Q'rl | ^t]- 
1=1 

By definition, 

_ > ^ I ^(^i) V ^oo) = exp{- / XJv} 

which implies that 

^P(r^+i -T^ >u\ a(r^) V G^) = -A,.+„exp{- "^^ A„cii;}. 
Hence the joint density function of {r^, r^, . . . , r*} given ^oo is: 

/ri,T2,...,T»|eoo(^i' • • • , ^i) = A(^At2 • • • Af. exp{— A„(it4 
Let {Z-' j^^Q be a sequence of continuous stochastic processes := {-^^}m>o such that 

nu 

:= / ZlX^q^dv, 
Jo 

where Z° = 1. Thus 

E[l{ri<s}PTi(lT^-^(lTi-'2 ■ ■ - Qt^ I ^oo] 

= •••/ qhqt2---(lu^iPtjT\T^,...,T^g^{ti,t2,...,ti)dtidt2...dti 

Jq^ Jo Jo Jo, 

= / A„p„exp{- / XydvjZl^^du. 
Jo Jo 

Let Hu = J^iZi Then by the monotone convergence theorem, 

Hu= [ iY.Zi + l)X„q„dv = I {H^ + l)X^q^dv. 
Jo Jo 

Note that a cadlag process has at most countably many discontinuities on a compact 
interval. Since {Xu}u>o and {pu}u>o are cadlag processes, Hu is continuous in m G 3?+ and 
differentiable except on a set with countable points which has measure zero. Consequently, 
except on an evanescent set, Hu satisfies the following Ordinary Differential Equation 
(ODE): 

— ^ = {Hu + l)A„g„ , 
du 
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and Hq = 0. 

Solving the ODE yields: 



= exp 1^ XvQvdv^ - 1- 

Therefore, 

oo 

P{r <s\gt) = Yl E[E[1{^.<s}Pt^<1t^-^<1t^-^ ...qri \ Goo] I 
i=l 

= YE[j^ Z*~^A„p„exp{-y^ Xydv}du\Qt] 
= E[ {1 + Hu)XuPuexp{- Xi,dv}du \ Qt] 

JQ JO 

= E[ A„p„exp{- / p^Xi,dv}du \ Gt] 
Jo ^ Jo 

= 1 - exp{- / puXudu}. 
Jo 



□ 



7.2 Proof of Lemma 1 

Proof: Note that the conditional expectation is on the set {r < s} and that the set 
{r > s} is an atom in Is. For any elements UJ ^ Qt and a) G Hg, due to the Markov 
property of the Poisson process, a) and {r > t} are independent conditional on u and 
{r > s}. Consequently, using a version of the Bayes' rule and Proposition 1, 

Eil{r>t} IGtynsV Is) = l{r>s}E{l{r>t} \ Gt ^ ns V Is) 

P{t > t\uj,uj) 



1 



{t>s} 



l{r>s} 



P(t > s \ uj,uj) 

P{t > t,UJ \ iiJ.T > s] 



P{UJ \ UJ,T > s) 
l{r>s}P{r >t\uj,T> S) 
P{T>t \ oo) 



l{r>s} 



l{r>s} 



P(r > s I cu) 
P(r > t I Gt) 



P{t>s\ Gt) 
^ exp{- Jq puXudu} 

^^^'^exp{- J^PuXudu} 

l{r>s} exp{- jlpuXudu}. 



7.3 Proof of Lemma 2 

Proof: For any s < t, using the results in Lemma 1, 

E[l{t>T} — 1{s>t} I J^s] = -£^[l{r>s} — l{r>t} I J^s] 

= l{r>s}(l - E[exp{- puXudu} I J^s])- 



□ 
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Note that conditional on Q^o and for x > s the density of default time is given by 
d 

— P(r < a; I r > s, ^oo) = PxK exp{- J p^X^du}. 

Hence, 

E[J^ PuXul{u<T}du I Ts] = 1{t>s}E[J^ PuXu'^luKTjdu I J^s] 

= l{r>s}E[j!. phXh exp{- PuXudu}U'!^ puXudu)dh 

+ exp{- PuXudu}{J^ PuXudu) \ J's] 
= l{r>s}(l - E[exp{- J^PuXudu} I J'sl)- 



□ 
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